Majorana fermions from Landau quantization in a 
superconductor— topological-insulator hybrid structure 
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We show that the interplay of cyclotron motion and Andreev reflection experienced by massless- 
Dirac-like charge carriers in topological-insulator surface states generates a Majorana-particle exci- 
tation. Based on an envelope-function description of the Dirac- Andreev edge states, we discuss the 
kinematic properties of the Majorana mode and find them to be possible to be tuned by changing 
the superconductor's chemical potential and/or the magnitude of the perpendicular magnetic field. 
Our proposal opens up new possibilities for studying Majorana fermions in a controllable setup. 
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Introduction. - In recent years, the possibility of ob- 
serving and exploiting Majorana quasiparticles in con- 
densed matter systems has attracted a great deal of at- 
tention [Il-Q • Originally these elusive fermionic particles 
were proposed as real solutions of the Dirac equation fsl . 
Potential condensed-matter systems hosting realizations 
of these atypical fermions include the v = b/2 quantum- 
Hall state 00], p- wave superconductors [8| such as stron- 
tium ruthena te |9 | , semiconductor-superconductor het- 
erostructures pi]Hl2| and the surface of topological insu- 
lators [1^ [3] among others. An experimental signature 
of localized Majorana excitations are zero-bias conduc- 
tance anomalies p"5l - [l8| that may have been measured 
recently Jl, J^l- 

A particularly interesting route towards realizing chiral 
Majorana modes involves hybrid structures of a topologi- 
cal insulator (TI) and an s-wave superconductor [2l|, [22] ■ 
The requirement of broken time-reversal symmetry and 
gapped excitation spectrum for the surface states in the 
TI could be fulfilled by proximity to a ferromagnetic insu- 
lator [i^l or Zeeman splitting due to a magnetic field [2l| . 
Here we pursue an alternative scenario where Landau 
quantization of the surface states' orbital motion in a uni- 
form perpendicular magnetic field is the origin of a gap 
and breaking of time-reversal symmetry. Besides avoid- 
ing materials-science challenges associated with fabrica- 
tion of hybrid structures involving three different kinds 
of materials, our setup also offers several new features 
enabling the manipulation of the Majorana excitation's 
properties. 

Figure [1] illustrates the proposed sample geometry. 
The massless-Dirac-like surface states of a bulk three- 
dimensional TI material (e.g., Bi2Se3 [l3|) occupy the 
xy plane. A planar contact with a superconductor in 
the half-plane a; < induces a pair potential for this 
part of the TI surface (henceforth called the S region), 
while a uniform perpendicular magnetic field B — Bz 
is present in the half-plane x > (the N region of the 
TI surface). We assume the magnetic field to be fully 
screened from the S part and neglect Zeeman splitting 



throughout 23]. The interplay of cyclotron motion of 
charge carriers in the N region with Andreev refiection 
from the interface with the S region results in the for- 
mation of chiral Dirac-Andreev edge states similar to 
the ones discussed previously for an S-graphene hybrid 
structure The quantum description of these edge 

channels reveals one of them to be associated with a chi- 
ral Majorana fermion mode with tunable velocity and 
guiding-center-dependent electric charge. 

Basic model for the S-TI hybrid structure. - To de- 
scribe single-particle excitations in the S-N heterostruc- 
ture made from TI surface states, we em ploy the Dirac- 
Bogohubov-de Gennes (DBdG) equation [ll [H] 



A*(r)ao 



A(r) ao 
THt,{y)T- 



*(r)=£*(r), (1) 



with the pair potential A(r) = e"^Ao8(— x) finite only 
in the S region, ctq = l2x2 denoting the two-dimensional 
identity matrix, and ^^0(1") = wf[p + eA(r)] • cr being 
the massless-Dirac Hamiltonian for the TI surface states. 
The position r = {x,y) and momentum p = —ih[dx,dy) 




FIG. 1. (Color online) (a) Schematics of the proposed sample 
layout, (b) The coordinate system, (c) The semiclassical 
picture of Dirac-Andreev edge states (solid and dashed lines 
represent electrons and holes respectively). 
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are restricted to the TI surface, cr is the vector of Pauh 
matrices acting in spin space. Furthermore, T denotes 
the time-reversal operator, — e the electron charge, and 
A the vector potential associated with the magnetic field 
B = V X A. The excitation energy e is measured rel- 
ative to the chemical potential /i of the superconduc- 
tor, with the absolute zero of energy set to be at the 
Dirac (i.e., neutrality) point of the TI surface states. The 
wave function Vl/ from Eq. ([1]) is a spinor in Dirac-Nambu 
space, which can be expressed explicitly in terms of spin- 
resolved amplitudes as ^' = itj,, — 



To describe the uniform perpendicular magnetic field 
in the N region, we adopt the Landau gauge A = Bxy. 
Then the momentum hq parallel to the interface (i.e., in 
y direction) is a good quantum number of the DBdG 
Hamiltonian, and a general eigenspinor is of the form 



*ng(r) 



e 2 



(2) 



Here is a Pauli matrix acting in Nambu space, uo 
the identity in spin space, and n enumerates the energy 
(Landau) levels for a fixed q. The spinors ^nq{x) are 
solutions of the one-dimensional (ID) DBdG equation 
n{q) $„g(x) = s <^nq{x), with 



H{q) = hvf 



' Tzi-i)dx + <Jy ' 



Tzq 



To -—- X kdyx) 



(3) 
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linear equations J2n' ^n"n' il) '^m ' = £ anq ' for the co- 
efficients, in the process defining the general EF matrix 



where the Tj are Pauli matrices acting in Nambu space, projecting onto specific basis states $„"o yields a set of 
and To the identity in Nambu space. A calculation sim- 
ilar to the one performed in Ref. [13 yields the spectrum 
of Landau-level (LL) eigenenergies Snq and explicit ex- 
pressions for ^nqix) in the N and S regions (see the Sup- 
plemental Material [1^, part A, for more details). When 
/i 7^ 0, a finite dispersion near q = is exhibited (see 
Fig. [2] that also shows the dispersionless Landau levels 
corresponding to /i = 0), which indicates the emergence 
of chiral Andreev edge channels [13, [1^ that are localized 
near the interface between the S and N regions. 

Emergence of the Majorana mode. - We now exam- 
ine in greater detail the states associated with the lin- 
early dispersing part of the n = level. A general 
symmetry property of the DBdG equation mandates 
that, for any eigenstate ^'„g with excitation energy e„g, 
its particle-hole conjugate S 4'„g in Nambu space [28| is 
also 8-11 eigenstcitc and has Gxcitation energy — £nq- 

This 

symmetry implies that the state with quantum num- 
bers n — and g = is its own particle-hole conju- 
gate, thus exhibiting the defining property of a Majo- 
rana fermion [l[, [2, [J. Explicit inspection indeed verifies 
the relation S5'oo(i") = — i5'oo(r) (see the Supplemental 
Material [25|], part B, for more details). 

While the Majorana state vE'oo(r) has a localized spatial 
profile in the direction perpendicular to the S-N junction 
(i.e., the x-direction), it is completely delocalized in the 
direction parallel to the S-N junction. To explore this 
Majorana mode further, we have developed a descrip- 
tion of single-particle excitations with q ^ based on 
an adaptation of the envelope-function (EF) theory com- 
monly applied in semiconductor physics [29|. A general 
eigenstate of the Hamiltonian ^{q) can be expressed in 
terms of the eigenstatcs at g = as 



(a) 2 

V2 

1 

#0 

-1 

-V2 

-2 

(b) 4 



2 

■5||o 

-Vl 
-2 



1 ' 1 ' 


' 1 ' 1 


11=1 








n=0 








n=-l 

1,1, 





-4 



-2 

4b 



1 ' 1 ' 


' 1 ' 1 










_n=6 

^1 

1,1, 









$„,(x)=^aif)<f„,o(x). (4) 

n' 

Inserting this expansion into the ID DBdG equation and 



FIG. 2. (Color online) Dispersion relation e^q of single- 
particle excitations at the S-N junction on a TPs surface 
and subject to a strong perpendicular magnetic field. Lan- 
dau levels with n = 0, ±1 are shown for Ao = ^hvY /Ib and 
fj, — OAIIvf/Ib [panel (a)], = LSSvp/Zb [panel (b)], and 
^ = (dotted lines). Ib = ^h/\eB\ is the magnetic length. 
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Hamiltonian 



dx [$„..o(a:)]^ {Hiq) ~ H{0)} ^n'oix) (5) 



Specializing to the ID DBdG Hamiltonian, we find 



n{q)-n{Q) = hvpgay 



(6) 



q=0 



which is independent of spatially varying quantities (the 
vector and pair potentials). Application of the Hellmann- 
Feynman theorem facilitates further calculation of 
the EF matrix Hamiltonian (see the Supplemental Mate- 
rial [23, part E, for more details). With mode-dependent 
velocities 



1 dSnq 

h dq 



q=0 



) 

(7) 



we finally obtain an expression where intra-LL and inter- 
LL couplings are neatly separated: 



d 



A „ 



/oo 
dx [$„"o(a;)]^ 
-00 



dq 



^n'q{x) 



(8a) 
(8b) 



q=0 



The EF approach is particularly useful to consider 
states with small q ^ as it provides a natural plat- 
form for a perturbative treatment. To lowest (i.e., linear) 
order in g, the dispersion of the n = Dirac-Andreev 
edge state is given by Egq = ^^o? + 0{q^), and the 
second-quantized particle annihilation operator for the 
state with n = and finite small q is found to be 



7? = Co 



q ^ AnQ (c„ 

n>0 



ci)+Oiq' 



(9) 



The Cj in Eq. ([9]) are single-particle operators associ- 
ated with the EF basis states $ jo (x) , satisfying the usual 
fermionic anticommutation relations, and cJ = Cq. The 
structure of the expression ^ embodies the relation 



(10) 



which fundamentally arises as a consequence of particle- 
hole symmetry in the DBdG equation. The identity ([TU]) 
is a hallmark of a Majorana excitation [22j, because it 
enables construction of a real-space particle operator [3l| 



7A(y)-A/'A/ dqe^'^y^g 



(11) 



that satisfies the condition 7a (y)^ — 7a(2/)- 

Majorana fermions are often also called 'half-fermions' 
because an ordinary (complex, Dirac) electron can be 
considered to be made up of two Majorana (real) fermion 



degrees of freedom with perfectly synchronized dynam- 
ics. Thus the emergence of independent Majorana exci- 
tations in an electronic system is often due to the frac- 
tionalization of an ordinary electron state, brought about 
by complex correlations in a material's electronic ground 
state. Naturally, a fractionalization scenario always gen- 
erates pairs of free Majorana excitations, and the phys- 
ical separation of the constituents of the pair has to be 
ensured to avoid their 'recombination' into an ordinary 
electron. The breaking of time-reversal symmetry due to 
the applied magnetic field eliminates the possibility for 
there to be a doubly occupied zero-energy state at a; = 
with guiding center q — 0. As particle-hole symmetry 
is preserved, such a state has to be a chiral Majorana 
excitation. There is no partner-complement along the 
NS-interface a; = to form a Dirac fermion, hence this 
excitation will be stable (e.g. against disorder) as long 
as a zero-energy state can exist. 

Electric charge of Dirac-Andreev edge states. As 
is the case for Bogoliubov quasip articles in a conven- 
tional superconductor [13, l32| - [35| . the Dirac-Andreev 
edge states are coherent superpositions of particles and 
holes in Nambu space. Consequently, their effective elec- 
tric charge generally has a non-quantized value 
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(12) 
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FIG. 3. (Color online) Efltective electric charge e*"'' of the 
single-particle excitations for the parameters used in Figs. [2^ 
and b, numerically calculated using the normal-side wavefunc- 
tions in Eq. (|12|) . For both sets of parameters, e'""' = as 
expected for a Majorana state. Notice also the linear variation 
of e*"") for small g</-\ 
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FIG. 4. (Color online) Dependence of the Majorana-mode ve- 
locity vo on system parameters. /3 = y/h\eB\vF / Ao measures 
the magnetic-field strength, /i is the superconductor's chemi- 
cal potential, and Ao the magnitude of the superconducting 
pair potential. 

that depends on the particular state. Figure |3] shows the 
numerically calculated charges e^"'-' of states from the 
bands n = 0, ±1 given in Fig. [51 

Using the EF theory discussed above, we can express 
the charge of modes with small q ^ as 

e("«) = Bnn + 2q^eJ2 An'n e„'n + 0{q^) , (13a) 
with the basis-function-related charge parameters 

en'n ^ y [^n'a{x)]^ CTq ® ^no{x) ■ (13b) 

By definition, e„„ = e*^"*^) is the effective charge of the EF 
basis state with quantum number n. The g-proportional 
coupling of basis states results in a linear q dependence 
of e„g for q < Ig^, in agreement with Fig. [3] 

For the n — Q (Majorana) mode, e*^°'~*^ = —^i^i) and, 
as expected, the state at q = has zero effective charge. 
A general superposition 5'o(r) = S|g|<A '^9 ^og(r) of 
Majorana-mode states has total electric charge 

= ^ (14) 

proportional to its mean energy (e) = X]|g|<A lo^gP^Oq- 
Thus the Majorana mode's energy relaxation will be asso- 
ciated with charge conversion, which is a general feature 
for bogolon excitations in superconductors 0, [s^ . Fur- 
thermore, as the superposition associated with the real- 
space Majorana-particle operator 7a(j/) involves states 



with ±q in equal measure, its total charge vanishes as 
is expected for a true Majorana fermion. An experi- 
mental verification of this property can possibly be done 
along similar lines as used in recent experiments on edge- 
magnetoplasmonic excitations in graphene 36] and open 
up possibilities for electrically manipulating this mode. 

Tailoring Majorana-mode properties. - An interesting 
feature of this realization of a Majorana mode is that 
its velocity can be tuned via the external magnetic field. 
This enables a potentially easier route towards tailoring 
this characteristic property than changing the magnetiza- 
tion of a ferroma gnet ic insulator in structures of the type 
discussed in Ref. [22- The mode- velocity dependence on 
the external magnetic field and the chemical potential 
is shown in Fig. H) Furthermore, the charge within the 
n = Landau level can be controlled via the chemical po- 
tential. This is evident by comparing the slopes of e'"'-* 
at (7 = in the two panels in Fig. |3l 

Conclusions. - We have demonstrated that a 
Majorana mode can be realized using the interplay 
of Landau quantization and Andreev refiection in a 
superconductor-topological-insulator hybrid structure. 
This mode is robust against disorder due to its chiral na- 
ture and associated uniqueness of the zero-energy state 
that is localized at the interface, and its velocity is tun- 
able by adjusting the strength of the magnetic field. We 
have developed a theory to describe the coupling of this 
mode to other Dirac-Andreev edge states and showed 
how this coupling results in an energy-dependent electric 
charge for Majorana-mode states. This feature should 
enable creation and manipulation of such excitations us- 
ing external gate voltages. 
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NCCR Quantum Science and Technology. 



Appendix A: DBdG eigenspinors in the N region 

To simplify the notation we measure all energies 
in units of hvp/ls and lengths in units of Is. The 
superconducting pair potential is in general, a com- 
plex quantity A = Age"^. To have a transparent 
picture of the relative phase between the electronic 
and the hole excitations we make the ansatz \I/(a;,?/) = 
e"iy{u^{x)e'*^^,ui{x)e''^^^,Viix)e-''^/^, -v^{x)e-''^/^f 
and obtain a one-dimensional DBdG equation whose 
solutions that decay for x — >■ oo in region x > have the 
form 
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/ -zCee-^(Ai + £)i?(^+e)2/2-i(a; + g)e^^/2 



Cte- 



(Al) 



with Ha{x) denoting the Hermite function [24]. The 
complete solution is then obtained by demanding the con- 
servation of particle current along the boundary. Similar 
to Ref. [H we find 

„ _ -iChili - £)H(^-e)2/2-i(-g) 

(A2) 

and the dispersion relation is given by the solutions of 



eiU+eiq)U^,{-q) + l) 
Ha2/2{q) 



(A3) 



The solutions e„ [q) of Eq. (jA3p can be labeled with a 
Landau level (LL) index n — 0,±1,±2, . . .. Figure 2(a) 
in the main text shows the zeroth LL {n = 0), along 
with the neighboring n — ±1 LLs when ji — QAhvp/ls 
and A = bhyp/lB- The dotted lines at 0, ±\/2 show 
the dispersionless LLs corresponding to /i = 0. The var- 
ious levels saturate at V2{hvF/lB)sgn{n)y/\~ri \ — ^ for 
q — > — oo and \/2{hvF /lB)sgn{n)y^\ n \ + /i for q ^ oo. 
This suggests that an interesting regime can be reached 
by increasing /i, so that ti = ±1 levels start contribut- 
ing to the low energy excitations of the system. When 



the chemical potential /i is finite (as measured from the 
charge- neutrality point of the free Dirac system), the LLs 
acquire a dispersion around q = that signals the exis- 
tence of Andreev edge excitations [H, [2^ . 



Appendix B: Verification of "too as a Majorana state 



For e = and g = we get Ch 
energy state can be expressed as 



iCe, and this zero 



*oo(a;) = Cee' 



I -i^iH^./2-i{x)e'^/^ \ 
H^y2{.x)e''^'^ 
^i^„2/2(a;)e-^*/2 



V Mi?, 



11' 



(Bl) 



The particle-hole conjugation operator is given by 
S = cTyTyK. where a and r are Pauli matrices act- 

ing on spin and particle-hole space respectively, and /C 
is complex conjugation. Straightforward verification im- 
plies SvI/qo = —i'^oo, and SS^I'oo = ^'oo- Hence the state 
^'oo is a Majorana fermion. 

Appendix C: DBdG eigenspinors in the S region 

In the superconducting region x < the solutions of 
Eq. (1) in the main text take the form [4^ 



^ix,y) = Ae^^ye'^^'-'e'^'-' 



( e-'^e 



V 



0/2 
-i(t)/2 



_g-i7giSg*0/2 
g-i</./2 



(CI) 



where the parameters 9,j,ko,K are defined by 

if £ < Ao 
if £ > Ao ' 




(C2) 



and 7 = sin ^ j^,ko = -^//i^"-^, k — sin(0). 

It should be noted that the solutions of the DBdG 



Hamiltonian given by Eq. (jCip . are approximate solu- 
tions valid in the regime /.t » Aq, e for all values of q. We 
are interested in the eigenstate at g = where Eq. (|C1[) 
represents an exact solution. 

The coefficients A,B,Ce, and Ch can be obtained by 
demanding the continuity of the wavefunction at x = 



[2J]. Due to translational invariance in the y direc- 
tion, the wavefunction matching can be done for each q 
separately. Thus we obtain (for q — and e — 0) 

^iA + iB = Ce{-i)fiHf,2/2-i{0), 



iA-iB = CeH^2/2{Q), 



A 



B = ChH^2/2iO), 



A-B = Chi-i)nH, 



(0). 



(C3) 



Once again we obtain Ch = iCe- The eigenspinor for the 
Majorana state for x > is given by Eq. (jBl[) . while for 
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a; < it is 



, (C4) 



where 



B = ^(-^iJ^./2_i(0) +^i/^./2(0)). (C5) 



representing the Majorana fermion in first-quantized 
form in the entire space. 



Appendix D: Calculation of the envelope- function 
Hamiltonian using Eq. (5) 

Using the explicit expressions for the eigenstates given 
by Eqs. (|Aip and (jCip we can write the EF matrix Hamil- 
tonian as 



Once again we can verify that the state is its own 
particle-hole conjugate. Thus we have obtained the state where 

I 



. ('^ + ^,i'o) +(^' + ^,i"o) -I 

Sn"n' = V^(£«'0)^^(en"0)[2 = TT {giSn'O, Sn"o) + 9{£n"0, £n'o)) 

-f G'(£„'o,en"o)2 = 7r(5(-£„.o,-e«"o) +.9(-£n"o,-£n'o))]- 



(D2) 



I 

The various functions N{r]), g{r],(), and G{r],() can be 
expressed using Gamma functions: 



Nin) 



r(i - (£i±!ili; 



r(i-i^) 



((m - ^)2r2(i - it^) + 4r2(i - (^)) (2r;r(i - i^^) + A^^i-^^i^ + - ^) 



V 



(.+.)^ (0)7^ + (a^ + ^(0)^/l ^ ^ ^ / g(.+c)^ (0)7^ + (m + C)g(.+c)^ ,(0)/ 

V 



(/i - ?7)g (^_^)2 (0) 



1-4 



where ip and F denote the Digamma function and the 
Gamma function respectively [4l| . 



r 



straightforward Dirac notation, we have 



Appendix E: Calculation of the envelope-function 
Hamiltonian using the Hellmann-Feynman theorem 



The r.h.s of Eq. (6) in the main text lends itself to the 
application of the Hellman-Feynman theorem. Using a 
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dn{q) 



I n'O) = (n"0 I 



g=0 



|-{H(<z)|n'g)}-H(<z)|-|A) 



g=0 



(n"0 I 
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n'O) + (n"0|{e„.o-H(0)} 
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FIG. 5. Velocity of the Majorana mode along y direction FIG. 6. Velocity of the Majorana mode along y direction as 

as a function p/Ao for different values of the magnetic field. a function of /3 for different values of /i/Ao. See the text for 

/3 = \Jh\eB\vF I defines a dimensionless measure of the further details, 
magnetic field. See the text for further details. 



The expression given in Eq. (jElcP neatly separates the 
intra-band and inter-band matrix elements in the first 
and second terms, respectively. Inserting into Eq. (5) in 
the main text and using the definition Eq. (7) in the main 
text yields Eqs. (8) in the main text. 

Naturally, the full matrix i/„"„'(g) is not tractable. 
However, we are usually only interested in the bands in 
the vicinity of the chemical potential, as the relevant dy- 
namics of the material will take place mostly in these. 
The coupling to more remote bands will affect the pa- 
rameters in any reduced-band model (e.g., the effective 
mass in conventional semiconductor band models, or the 
velocity of the edge modes in our current system of in- 
terest). The systematic way to derive a reduced-band 
model is Lowdin partitioning (see, e.g.. Appendix B in 
Ref. 43] for a pedagogical introduction). 



Appendix F: Majorana mode coupled to nearest 
bands as an illustration of envelope-function theory 

We will now consider the coupling of the Majorana 
mode (n = 0) to the nearest LLs (with n — ±1). The EF 
matrix Hamiltonian for this situation is given by (we use 
the notation 1 = — 1) 

(hvoq SiAoiq -eiAoiq \ 
-eiAiog ei + hviq -2eiAm , (Fl) 
eiAiQq 2eiAiiq ~ei + hviq J 



where Si = £io = — ^lo- Straightforward diagonal- 
ization yields eigenvalues £±1, = ±£i -I- hviq + 0[q^), 
eoq — hvoq + 0{q^) and associated eigenstates Xnq = 

/ 1 \ ( Aoiq\ ( Aoiq\ 

XOq= Aloq ,Xl9= 1 >Xl9= All , 

\AioqJ \AuqJ \ 1 J 

(F2) 

also with accuracy of 0{q^). Thus, up to linear order 
in q, the inter-band coupling does not affect the linear 
dispersion of each mode, but an admixture between the 
different modes is present in the eigenstates. 

We now consider the second-quantized representation 
of the Dirac-Andreev edge states. Based on the above 
discussion, we can write the annihilation operator for the 
state with wave number q in the n = band as 

lq = co + q {AiQ ci -I- Ci) . (F3a) 

Here c„ are the second-quantized operators associated 
with the q = states in band n, i.e., our basis states. By 
the symmetry inherent to the Dirac-BdG equations, we 
have Co = cj and Aiq cj = —Aiq c\ and thus find 

-fq = Co + q Aio{ci - c\) . (F3b) 
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Appendix G: Tuning the system properties via the 
external magnetic field 



It should be noted that all the energy scales and mo- 
mentum scales in Eq. (jA3l) are defined using the magnetic 
length thus removing any explicit dependence upon the 
external magnetic field. If we are interested in tuning 
the properties of our system by changing the external 
magnetic field, it is useful to define another set of di- 
mensionless units where we can measure the strength of 
the magnetic field, this is achieved by measuring all the 
energies in units of the superconducting gap. We define 

= ^ ^ Ao?' ^ ^ a'^^ where /? is a dimensionless 



measure of the magnetic field (/3 
Eq. (|A3[I can be rewritten as 



y/h\eB\vF y Using this 



, (Gl) 



where the function fa{q) is the same as defined earlier. 
Equation (jGl[) can be solved for e. The numerical results 
for the velocity of the Majorana mode as a function of 
/? and /i/Ao is shown in Fig. 4 in the main text. In this 
supplemental material Fig.[S]shows vq/vf as a function of 
/i/Ao for different values of /3 and similarly Fig. [6] shows 
vq/vf as a function of /3 for different values of fJ-/ ^.q just 
to illustrate the oscillating behavior of vq/vf as seen in 
Fig. 4 in the main text. 
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